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Abstract 

We study classical dynamics of a probe Dp-brane moving in a background sourced by 
a stack of Dp-branes. In this context the physics is similar to that of the effective action 
for open-string tachyon condensation, but with a power-law runaway potential. We 
show that small inhomogeneous ripples of the probe brane embedding grow with time, 
leading to folding of the brane as it moves. We give a full nonlinear analytical treatment 
of inhomogeneous brane dynamics, suitable for the Dirac-Born-Infeld + Wess-Zumino 
theory with arbitrary runaway potential, in the case where the source branes are BPS. 
In the near-horizon geometry, the inhomogeneous brane motion has a dual description 
in terms of free streaming of massive relativistic test particles originating from the 
initial hypersurface of the probe brane. We discuss limitations of the effective action 
description around loci of self-crossing of the probe brane (caustics). We also discuss 
the effect of brane folding in application to the theory of cosmological fluctuations in 
string theory inflation. 
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1 Introduction 



In this paper we investigate generic inhomogeneous solutions describing a probe Dp- 
brane or anti-Dp-brane falling in the near-horizon D-brane background geometry. A 
particularly interesting case occurs for D3-branes, where the background involves anti 
de Sitter space. The background geometry is produced by the stack of N source Dp- 
branes. We will show that dynamics of brane motion can be reduced to that of the 
Dirac-Born-Infeld (DBI) type, with the action 



Here T(x fl ) is a scalar field (dimensionless: in units of a') and V(T) is its runaway 
potential with no minima. The action should be understood in the truncated approx- 
imation. While first derivatives can be arbitrarily large, the action is valid only in the 
regime where second and higher derivatives are small. The action (^Q) is identical to 
the effective field theory action of the open string theory tachyon describing unstable 
non-BPS D-branes as found in e.g. 1 HI Ell- The similarity between the DBI effective 
action for the tachyon and that of a probe D-brane near NS5-Branes, for homogeneous 
solutions, was recently studied by Kutasov 4 J. Some other recent works on similar 
topics (again for the homogeneous case) include El El • 

Previous works |10|lll| on the generic inhomogeneous solution in the theory Q were 
focused on the open string theory tachyon effective action, with runaway potential 
of exponential asymptotic form: V{T) ~ exp(— T). A generic analytic solution was 
found jllj which, in particular, shows the growth of T inhomogeneities resulting in 
tachyon energy fragmentation and folding of T. In our paper here, we extend the 
method of JJ| to treat generic inhomogeneities for a test D-brane or anti-D-brane 
freely falling in a background D-brane geometry. Folding of the field T in this context 
acquires a transparent interpretation, as the folding of the brane embedded in the 
higher dimensional spacetime background. Folding of an inhomogeneous probe D-brane 
should also transpire for the background geometry generated by source NS5-Branes. 

The setup of moving a brane in a background spacetime generated by many branes 
is often an ingredient of string theory cosmological construction including models of 
inflation, as in e.g. |131 114j . Usually generation of inflationary perturbations in these 
models are treated in the framework of four-dimensional effective theory. The effect of 
the brane folding which we find here is based on the full treatment of the probe brane 
dynamics. Thus it should be incorporated in the complete treatment of fluctuations in 
cosmological applications. 

Before proceeding with the story of an inhomogeneous probe D-brane in the back- 
ground geometry, let us recall the story of the inhomogeneous rolling tachyon. Study 
of the rolling homogeneous tachyon in classical string field theory has shown that, at 
late times, the system is described by a pressureless fluid "tachyon matter", with no 
perturbative open string modes at the bottom of the potential. Turning on g s would 
lead to subsequent decay of the unstable brane, mainly to massive closed string modes 
as in [15] . According to the proposal of |T§1, this process has a dual description purely 
in terms of the rolling open string tachyon. 

It is difficult to find exact solvable CFTs describing open string tachyon dynamics 
for a generic inhomogeneous tachyon profile. Some progress in this regard was made 

We do not attempt to be exhaustive referencing the string field theory literature here. 
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concerning evolution of the string theory tachyon with a plane wave tachyon profile • 
In this case, the tachyon decays into equidistant plane-parallel singular hypersurfaces 
of codimension one, which were interpreted as kinks. This inhomogeneous profile is 
atypical, though, in the sense that fragmentation between kinks does not occur. In 
the general case, based on cosmological intuition, we expect both types of structures: 
weblike fragmentation and topological defects. 

Quite apart from progress in the CFT approach, the action Q has in fact proven 
to capture many striking features of tachyon condensation. The relatively simpler 
formulation of tachyon dynamics, in terms of the effective action (^Q), has therefore 
triggered significant interest in investigation of the field theory of the tachyon, the 
possible role of tachyons in cosmology, and so forth. Indeed, the end point of string 
theory brane inflation is annihilation of D-branes and anti-D-branes, which leads to 
the formation and subsequent fragmentation of a tachyon condensate (101 lllj . So 
the potential role of the tachyon in cosmology cannot be understood without first 
understanding its fragmentation. 

In |10| . it was found that the evolution of the tachyon field T(t, x) can be viewed as 
a mapping T(to,xo) — ► T(t,x), that evolves to become multi- valued, with singularities 
at caustics generated. This seemed to be a generic behavior for runaway potentials. 
More recently, a study considered the generic inhomogeneous tachyon field, in the 
region where it rolls down one side of the potential. Formation of sharp features in 
the tachyon energy density, due to fragmentation, was observed. The tachyon energy 
density pattern is reminiscent of the illumination pattern at the bottom of a swimming 
pool, or the web-like large scale structure of the universe. The similarity is not coin- 
cidental: the underlying mathematics has common features in all three cases. In the 
free streaming approximation of 12] this is just focusing of particle trajectories corre- 
sponding to higher density concentrations and, further, to the formation of caustics at 
the loci where trajectories cross. These features, which are related to the convergence 
of characteristics of the field T, are distinguished from topological defects. The full 
picture must incorporate both effects: formation of kinks and tachyon fragmentation in 
the space between them. In the context of brane inflation ending with annihilation of 
D-branes and anti-D-branes, the tachyon is a complex field and strings will be created. 
The web-like structure of fragmentation which will be superposed with the network of 
strings. 

In this paper, we look at another situation, concerning the radial motion of a probe 
brane in particular string theory backgrounds. Interestingly, moving probe D-branes 
in the background of RR-charged or solitonic branes of various dimensions gives rise 
to an effective probe dynamics in which the DBI part, after some field redefinition, has 
the same structure as (0) - with a runaway potential of power law type, i.e. 

V(T) = cT- a (2) 

for some a > 2 . For example, in the recent work 0], dynamics of a BPS test D-brane 
in the vicinity of stack of NS5-branes was considered, and for radial motion of the 
probe an action of type Q arises, with a 'tachyon' potential falling off exponentially 
for large value of the 'tachyon' field. In general, the DBI part of the probe action 
for BPS Dp-branes probing a stack of other branes in the near-horizon region of the 

2 Clearly, this behavior differs from that of the usual open string tachyon potential, which is exponentially 
suppressed for large values of tachyon field. 
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background has dynamics of type (0), but with a runaway potential (J2J) with index 
a = 2(7 — p)/(5 — p), which is positive, and hence sensible, for < p < 5. 

Clearly, in these cases, there is a spacetime interpretation of the 'condensation' 
of the radial 'tachyon' mode, in terms of motion of the probe brane in a background 
geometry. Analysis of such may open a window to a better understanding of the full 
open string tachyon condensation itself [3]. 

In cosmological applications, the radial distance r to the probe brane plays the role 
of the inflaton field ^21- in a particular cosmological setup, cosmological scenarios 
based on the motion of a test anti-D3-brane in the throat region of a large number of 
D3-branes have recently been studied |13[ I14j . Dynamics of the radial mode 3 is given 
by an action of the above type (0) with a power-law runaway potential with a = 4. 

Inhomogeneous fluctuations of this field during inflation generate primordial cos- 
mological fluctuations. Therefore it is important to know the time evolution for not 
only homogeneous, plane parallel brane motion in the anti de Sitter or more general 
background geometry, but generic inhomogeneous branes. 

Understanding of the generic solution of the DBI type effective action may give 
insight to the structure of the theory. 

The plan of our paper in the following. In Section [21 we introduce the supergravity 
description of the probe brane freely falling in the background geometry. We will 
show that the brane dynamics is reduced to the compact form of the DBI plus WZ 
effective action for a single scalar T with the power law potential. In Section |H1 we 
make the first approximation towards the inhomogeneous solution of the non-linear 
equation of motion for T, the free-streaming approximation. Yet, the free-streaming 
approximation is not enough to calculate the stress-energy tensor of the brane. In 
Section0]we construct the full analytic solution in terms of asymptotic series expansion 
for T. Most importantly, the next-to-leading term allows to calculate the stress-energy 
of the brane. The answer is transparent and allows a conjecture of a duality between a 
probe brane freely falling into the near-horizon background geometry (such as the AdS 
throat of D3-branes) and a collection of massive relativistic particles, see Section |21 In 
this final section we also discuss various aspects of the physics of the folding branes 
and its cosmological applications. 

2 Probe brane in background brane geometry 

The physics of a probe Dp-brane or anti-Dp-brane moving in the background of N 
source Dp-branes is given by the usual string-frame DBI plus Wess-Zumino action, 

Sprobe = Sdbi + S wz = -Tp j dP +1 ae~^^- det(P(<7 a/3)) =p r p f P(C (p+1) ) (3) 

where =F denotes - for probe anti-branes or + for branes, and P is for pullback. 

For a (non-extremal) Dp-brane spacetime in string frame, in a standard coordinate 
system, 

e <£ = H (3- P )/4 

ds 2 = H- 1 l 2 (-Kdt 2 + dx\^+H 1 l 2 (dr 2 /K + r 2 dnl_ p ) 

C( p+ i) = C l (1 - H' 1 ) dt A dx\ • • • A dx p (4) 
3 Apart from its coupling to the 4D gravity, that is. 
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Figure 1: Sketch of probe brane moving in the background geometry. 



The harmonic function is given by 
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It follows easily that 
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where for convenience we use an auxiliary metric (guv) = 1, ... 1, 1/K, r 2 . . . r 2 ), 

which collapses to the flat metric in the BPS limit. In static gauge, where a a = x a for 
a = 0, . . . ,p, this expression simplifies in the usual way to give 
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(static) 
probe 
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1 + Hg ij d a j 1 dpx?jr ifi ± C 



(9) 



A sketch of the probe brane in the background sourced by a stack of N Dp-branes is 
shown in Figure ^ 
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Now let us restrict to configurations in which the transverse motion is purely radial. 
In this case, the physics simplifies to that of just one field describing the motion of 
the test (probe) brane. This can be seen by performing a field redefinition to give a 
'tachyon' 

T ^I dH ^ H - 1) - <s - pW -' ) - (10> 

Physically, this expression (|1U|). in combination with the expression © for H(r), gives 
the explicit connection of the radial coordinate r and the field T. In other words, 
T has clear geometrical interpretation in terms of the radial position of the brane 
embedding as r = r(t, x\\). Mathematically, in general, the relation (jl()j) gives T in 
terms of hypergeometric functions of H. 
The field redefinition yields 



s (static,radial) = _ r J jp+l^y^ 



l + g^d a TdBT±( 



(11) 



where the indices run only over a, (5 = 0, . . . ,p and the only nontrivial component of g 
is the time component. In this probe action (jllj) for radial-only motion, the tachyon 
potential is given implicitly by 

v{T ^mry <12> 

To obtain the explicit form of V(T), it is necessary to invert the T{H) expression (|1()|) 
to get H{T). For Dp-brane geometries where the probe is close to the stack of source 
branes (e.g. in the decoupling limit in the string theory context) 

(13) 

the above inversion can be done relatively easily. We assume this holds in the following. 
For the general case for source branes (BPS or not), with the definition 

we obtain the exact expression for T{r) in terms of hypergeometric functions: 

T(r) = 2fi f I -; i + -; i - K(t)) (15) 

Mathematically, this expression is monotonic in r, and is convergent for K(r) £ [0, 1] 
(including at the endpoints) because a > 0. 

When the probe is far from the horizon, in the sense that K(r) ~ 1~, the series 
expansion of the hypergeometric function can be employed; this is obviously also valid 
when there is no horizon (BPS source branes). We then get 
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Hence, large-T corresponds to small coordinate distance from the source branes. Since 
we are considering the large- H limit (|13[). this is always our regime of interest. The 
potential for the tachyon 



2R 
(5~p) 



V(T)\ hps = (j^xT) (17) 



exhibits power-law runaway behavior. 

In the non-BPS case, in order to obtain the near-horizon K{r) ~ + expansion of 
(|15|) . we use a standard identity for hypergeometric functions to get 

.fi(j-;i + -;i-K) = ffi±^f 2 F l( i/2,i/ a ;i/2^) 

\2 a a J 1(1/2 + 1/0) 

~ 2 r(t/a) ^2^1 (1/2 + 1/a, 1; 3/2; if) (18) 

and then, again, employ the series expansion. The upshot is that, in the near-horizon 
limit, the entire expression is finite all the way to the horizon. The second term in 
(|18jl is subleading at small- K, while the first term has a coefficient as a ratio of two 
finite T-functions. Amusingly, then, the non-BPS expressions for T{H), and hence 
V(T), end up being identical to the BPS case, up to a finite numerical multiplicative 
constant. This system therefore exhibits the same qualitative physics, as far as the 
runaway tachyon potential is concerned. The physically relevant qualitative difference 
for non-BPS branes is, of course, that the auxiliary metric on the brane worldvolume, 
iSap) = diag(— K (T), 1, • • • ,1) features a non-trivial time slowdown near the horizon. 
In the following, we will choose to stick to the BPS case, because it can be solved 
exactly in an analytic asymptotic series expansion. 



3 Free-streaming approximation 

The equation of motion for the tachyon field follows from the action (|11|) 



UT ~ i + d a Td«T ^ ryr ~ ir V + CVl + d*Td a T) = . (19) 

This equation is an example of a non-linear, partial differential equation which, as we 
will show, admits a relatively simple, general, inhomogeneous solution. 
The energy density of the tachyon field p = Tq is 

V(T) 



y/1 + d^Td^T 



T + V(T) y 7 ! + d^Td^T, (20) 



where V = V(T) is given above in (|17|). As has been observed in ^U] by solving (|19|) 
numerically, if we define an operator P(T) = l + d^Td^T the field T rapidly approaches 
to a regime in which P(T) ~ 0. 

Importantly, this means that in the same regime, the Ramond-Ramond coupling is 
essentially irrelevant physically. In particular, neither the ± sign (whether the probe 
is a brane or an anti-brane) nor the size of £ (the non-extremality parameter) matters! 
This conclusion is somewhat surprising, naively, because one would think that the 
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Ramond-Ramond fields would be very important near the source branes. The point is 
that the gauge fields become essentially irrelevant for this test-brane falling close in to 
the source branes, because of a strong kinetic suppression. We will nonetheless keep 
track of the R-R contribution to the theory on the probe. 

A summary of the free-streaming tool of JOJ goes as follows. In the leading ap- 
proximation, 

T(X") « S(X") , (21) 

where S satisfies the equation 

S 2 - (v x Sy = l, (22) 

where we use t = X° , dot is dt and the spatial derivatives are with respect to the 
p spatial coordinates x = on the brane. This equation is the Hamilton-Jacobi 
equation for the evolution of the wave front function of free steaming massive relativistic 
particles. In this particle description, at some initial time to we can label the position 
of each particle with a vector q\ equivalently, we can say that q parameterizes the 
different particles. The initial (covariant) velocity of the particle is given by O^Sq. If 
we further define the proper time r along each particle's trajectory, we can switch from 
coordinates (t,x) to (r(t,x),q(t,x)) and obtain an exact parametric solution to (|22|) 

m 

x = q-VfS T 

t = ^l + |V g -5 | 2 r 

S = S + T. (23) 

Notice that the coordinates (r, q) are the proper coordinates of the probe brane. The 
interpretation of the solution (|23|) is very simple and intuitive. It tells us that the 
field S propagates along the trajectories of the massive relativistic particles, growing 
linearly with proper time. The slope of each characteristic depends only on the initial 
gradients of So on that characteristic. 

In the free streaming approximation, the denominator of the first term in the 
tachyon energy density expression (|2U|) is zero, because the leading term in T 2 is unity. 
Given this fact, it is clear that in order to compute the energy density near the asymp- 
totic runaway regime, it is necessary to go beyond the free streaming approximation 
(|22jl. In the next section, we compute the energy density of the tachyon by finding 
corrections beyond the free streaming approximation by using a working ansatz in an 
asymptotic series expansion. 

4 The full solution 

We are interested in the energy density: something the free-streaming approximation 
cannot give us. So let us begin here by describing the energy density qualitatively. 
Looking at equation ()2()D we see that the potential pieces are growing small, as an 
inverse power-law, as are the arguments of the square roots. The second term in (|20(l 
will thus rapidly become irrelevant. What we need consider is only the competition 
in the first term in (|2U|) between the small potential in the numerator and the kinetic 
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square-root term in the denominator, all evaluated near the runaway asymptotic re- 
gion of the potential. In a previous work jllj . the asymptotic series solution for the 
exponentially decaying potentials was proposed to be T(x M ) « S + XmLo fne~^ n+1 ^ S ■ 
Here we generalize this result by conjecturing the following asymptotic expansion for 
the field T 4 

oo 

T(x") -S + Y, f n S- {n+1)l (24) 

n=0 

In the next to the leading approximation, we have 

T(X») + fid , (25) 
where G\ is of the following form: 

d = S~ 7 , (26) 

where 7 > and 

» = ^1 (27) 

Plugging the expansion (|25|) into equation (|19() and keeping only terms linear in /1, 
gives the following equation for the function f\ 

- sd»sd u s{d^d v h) + (2sas - 2 7 + 2a) d»s(d li f 1 )+ 

[(2a 7 -7(7 + l))5- 1 +2 7 D5]/ 1 =0, (28) 

where OS = d v d u S = —S + V^ 2 ^. This equation can be dramatically simplified by 
changing from it, x) coordinates to (r, q) coordinates, as defined by the characteristics 
of S in (|23j) . In these coordinates, /i has no spatial derivatives and equation (|28|) 
reduces to 

Sf 1>TT + (2SDS -7 + 1) h,r ~ ^US h = 0, (29) 

where using (|27j) we have got rid of the term proportional to S . Note that OS can be 
calculated either with respect to (t, x) or (r, q) coordinates. We can further simplify 
this equation by introducing a new variable y 

y = ih- Sh,r (30) 

The second order differential equation (J29|) can be rewritten in terms of y satisfying a 
first order differential equation 

y >T + 2nSy = 0, (31) 

which can be immediately integrated to give y(r, q) = yo(q) exp (—2 f T dr' OS) . From 
this and (fHOjl 

fi(r,q) = Mq)Sl fdi'Stf, fr^e' 2 ^ dT " ns . (32) 

(The integration constant that one gains after integrating (|30|) is irrelevant as it can 
be absorbed into S.) 



4 One can try asymptotic series with respect to other of S. However the final answer shall be independent 
of the choice. 
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To proceed further we need to estimate 1)32(1 . It can be shown that exp(— 2 j T dr'nS) 
is a polynomial of order 2p with respect to r such that f\G\ is small in the asymptotic 
region. We start by referring back to the free-streaming parametric solution ()23j) . 
Using the fact that T ~ S + fiG\, this shows the validity of the asymptotic series 
expansion (|25|) for runaway power-law potentials. We discuss second-order corrections 
in the Appendix A; these are also shown to preserve the validity of the asymptotic series 
expansion. Comments on arbitrary runaway potentials can be found in the Appendix 
B. 

To compute the energy density, we need more information. Also, proper coordinates 
(t, q) are not so easily interpreted in the spacetime picture, where the brane position 
in static gauge is given by X a . Therefore, let us compute the Jacobian that transforms 
( T ; 0) ~* (£> To see its role, let us use the following trick. Plugging the expansion |2f 



into the energy density (|2U|) we find p » 1/ (7/1 + Sd^Sd^fi). The denominator 
here is none other than \/2y. Therefore, 



2y 



Po (q) exp dr'ns) . (33) 



This is precisely the expression for the energy density of free streaming relativistic 
particles which obey the relativistic continuity equation d^S (pd^S) = 0. In fact, 1)33(1 
is the solution of this continuity equation in the coordinates (r,q). However, there is 
another form of the energy density, which is equivalent to ()33j) 



P = Po(o)/ 



dx 



dq 



(34) 



where the denominator is the Jacobian of the x — > q transformation. Note that we can 
compare expressions ()3*3*)l and for p to find a precise expression for \3S and then, 
via ()32)) . for f\. Indeed, from conservation of the energy density in a differential volume 
we have d p qpo(q) = d p xp, which leads to the formula ()34j) . It is now straightforward 
to calculate the Jacobian from the formulas ()23j) . In p dimensions it is a polynomial in 
r of order p. This can be seen by using x = q — rV q So. We can define the symmetric 
tensor D£ = V^V^So; then the Jacobian can be written in terms of the eigenvalues of 
D, = 1. . .p}, as the product (1 — Air) ... (1 — A p ). Using ()34j) for p, it therefore 

becomes clear that there is a critical proper time where the energy density blows up. 
In detail, at leading order in the asymptotic series the energy density is given by 

p 1 

f=i 0- - a<t) 

Corrections to this expression are suppressed by powers of S. 

A similar behavior in the context of exponential potentials was observed in 
We now discuss this phenomenon in our new context of power-law runaway potentials, 
relevant in string theory for D-brane motion, in the next section. 



5 Interpretation and Discussion 

Consider the probe brane or anti-brane in the background geometry. Suppose the 
generic embedding is not exactly located in the hyperplane Xu, but rather given by 
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Figure 2: Folding of the brane segment from initial profile 1 through the instance 2 of the 
first caustic formation to typical multiple valued configuration 3. 



the distance function r = r(xu). Then the following picture emerges of how the brane 
is moving in the background geometry. First, very quickly the brane hypersurface 
becomes moving in the regime of free streaming, when each point of the brane is 
ballisticaHy propagating. By very quickly we mean that the test brane is in the vicinity 
of the near-horizon throat. Simple formulas Q23|) . resembling the Huygens principle of 
geometric optics, are describing the evolution of the brane hypersurface T; in other 
words, its radial position r. Geometrical rules and pictorial illustrations of the free 
streaming approximation are given in detail in the previous papers |1U1 lllj . 

The most important feature of brane evolution we see here is unbounded bending 
of its hypersurface, in other words the growth of its inhomogeneities. A simple rule 
stands, that the brane regions around local maxima of T are getting flatter while 
regions around minima are sharpening. This illustrated by the Figure |2] which shows 
edge-on profile of the folding brane segment around the minima of T with the most 
dramatic reconstruction. At some moment the brane profile acquires a discontinuity 
where its second derivative blows up. This corresponds to caustic formation. Formally, 
if we allow the free streaming approximation to work further, the brane hypersurface 
is self-crossing, and T goes multi valued as shown on the figure. 

One may think that the multi- valuedness of the field T is related to the failure of 
the parametrization of the world-volume of the probe Dp-brane in static gauge. We 
argue that brane self-crossing is real, and accompanied by local blowup of the stress- 
energy tensor. Indeed, we went beyond the free-streaming approximation to describe 
evolution of T. While next-to-leading order terms in the series T ~ [S + 0(1/S n )] have 
essentially negligible correction to the free-streaming propagation of T, they are fully 
responsible for the stress-energy tensor. However, the very simple and transparent 
expression for stress-energy tensor (|33l34j) shows that at a certain critical time and 
location the energy density becomes singular. The formula (|33I34|) is given in terms 
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of proper time and proper brane coordinates, and does not relate to the static gauge. 
We conclude that formally the effective DBI+WZ theory for generic inhomogeneous 
solutions for radial probe-brane motion leads to multivalued solutions, with caustics 
where the stress-energy tensor turns singular. 

Certainly, the effective field theory description of the brane dynamics is limited by 
finite field gradients and stress energy components. Therefore, we have to discount 
the formal results beyond the points where the field gradients are large (say, larger 
than unity in string units, which we used here). What happens with the probe brane 
around a region about to develop a caustic? Here, we only partly address the issue 
and suggest some directions for future thoughts. One possibility is related to the 
classical physics of extended objects with sharp kink-like features. Sharp restructuring 
of the brane hypersurface may lead to significant radiation similar to gravitational wave 
flashes from cosmic string kinks. This would result in local dissipation of the stress- 
energy into classical radiation. Another possibility is better couched in terms of string 
theoretic physics. Suppose the multi- valued configuration as shown in Figure [21 begins 
to form. The horizontal segment of brane in configuration 3 has orientation opposite 
to that of the rest of the brane. That is to say, this segment now effectively describes 
a piece of anti-brane. Brane-antibrane segments in close proximity are unstable due to 
the open-string tachyon - corresponding to fundamental string stretched between the 
segments - which results in their annihilation. 

Next, we discuss implication of our effect for cosmological applications. 

The inter-brane moduli field r is often exploited as an inflaton in the string theory 
constructions of the brane inflation. We mention one of the realization of this idea 
which is close to the setup of paper, namely the inflationary model of ^3]. The high 
dimensional construction includes a D3 brane freely falling in the background of five- 
dimensional anti de Sitter geometry generated by a large stack of source D3-branes. The 
radius r plays the role of the moduli field in the effective four-dimensional description. 
Under certain conditions, again, in the 4d effective description of gravity plus that 
moduli field inflationary regime can be realized. 

One of the most important prediction of inflation is generation of scalar metric 
perturbations, usually associated with quantum fluctuations of the inflaton. Therefore 
we are interested in spatial inhomogeneities of the moduli field. In the context of our 
setup, these correspond to the spatial inhomogeneities of the brane embedding. 

As we learned, even at the classical level brane inhomogeneities are unstable and 
grow. It is therefore interesting to consider combining probe-brane dynamics with 
inflation in 3 + 1 dimensions. At the level when 4-dimensional gravity is introduced 
phenomenologically (corresponding to adding an Einstein-Hilbert term in the action 
by hand), in the DBI effective brane action we just have to introduce four-dimensional 
(quasi-) de Sitter geometry. In the equation of motion for T, this will just lead to the 
friction term 3HT and redshifting of all spatial gradients — ► ^d^. We can conjecture 
that apparently folding of the brane will be suppressed at scales of the de Sitter radius 
1/H. However, it can persist at larger scales. A more quantitative description of this 
is left for the future. 
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Appendix A 



We can use a trick to derive the second-order terms in the perturbation about free 
streaming. Namely, we write 



fiG 1 ^G 1 (f 1 + f 2 G 2 ) 



(36) 



where G 2 is a power-law function G 2 cc S~ x , whose index A to be determined. 

The above trick contains nontrivial information. Indeed, for any n-th order pertur- 
bation f n we will need to use the equation of motion for the lower-order perturbation(s), 
which in general yield an equation for f n which is not identical to the equation for f\. 
Let us see this explicitly. The equations give 



{d li d v f 2 ){d^WS)S - (d^WS) [(2a - 2 7 - 2A) + 2SDS] 

A(A + 7 ) 



h 



2( 7 + A)QS- 



S 2 



= 



2( 7 + A)OS- 



Changing coordinates to (r, q) gives 

5/ 2)TT + [25D5+(l-7-2A)]/2, 
Defining 

V = Sf 2 , T G 2 + ( 7 + X)f 2 G 2 
we find that again the equation becomes integrable: 



A(A + 7 ) 
S 2 



(37) 

(38) 
(39) 



S(t>) 




exp - 



S(t") 



(40) 



At long-time, S ~ r, and this equation can be integrated easily. It is then clear that 
f 2 G 2 is of the same order as f\. In other words, the extra piece just renormalizes f\. 
This fact holds true for any choice of A. 

We can obtain higher-order equations in perturbation theory by repeating our trick 
iteratively: substituting f n -i — ► G n -i(f n -\ + f n G n ) and using lower-order equations of 
motion. We conjecture that at general order additional perturbations either correspond 
to a reparametrization of S, or a renormalization of lower-order terms, or pieces which 
are down by additional powers of S in an asymptotic expansion. 
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Appendix B 



It is interesting to ask if the analysis at first order beyond the free-streaming approxi- 
mation can be done for arbitrary runaway tachyon potential. Indeed, we show how to 
do this here. 

We start by writing 

T(x»)ttS + e(T,q) (41) 

Next, we plug this expression into the general nonlinear equation of motion (|19|) for 
the tachyon. 

In an asymptotic series expansion, e should be much smaller than S, this will prove 
to be a self-consistent approximation. Keeping only lowest-order terms in e, and using 
the equation of motion for S, we find a surprise: the equation for e involves only 
derivative terms. 

The next simplification occurs upon changing to proper coordinates (r, q) . The 
resulting equation of motion for e becomes 



+ e,r 



-205 + 2 



T=S- 







(42) 



which can be easily integrated (twice) to give 



((T.rf) - (,(<]) I ,/r"exp ( 2 / dr'^jexp 



dr'DS 



(43) 



In the long-time limit, or when initial inhomogeneities are small (which is valid 
here) we can do the integral involving the potential to give 



<r,q) = eM / dr" [V (S(r"))] 2 exp -2 / dr'DS 



(44) 



As is clear from the above equation, in the long-time limit, the suppression of the next- 
to-leading term is controlled by the square of the tachyon potential, as was suggested 
in section 4. For runaway potentials, this is always small in the regime of validity of 
the asymptotic series expansion. In particular, note that our results here agree with 
the case of the power-law worked out in the body of the paper. 
The energy-momentum tensor can then be written as 



V(S) 



(45) 



where 



A J2d@ Sdpe + (higher order terms) 



(46) 



We have d^Sdpe ~ F 2 exp(— 2 f T dr'DS). Substituting this into equation (|43)) we 
obtain equations (j3*3*|) . (|3*1|) . Thus we find out that equations (j3"3*|) and (|3*1|) for the 

energy density are valid even for arbitrary runaway potentials. This quantity vA 
controls the degree of kinetic suppression, of e.g. the Ramond-Ramond terms in the 
D-brane probe action. 



'In proper coordinates, of course, d"Sdp£ — —2e T . 
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